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ABSTRACT 


Nonequilibrium  performance  losses  are  interpreted  in  terms  of 
radical  concentrations  and  extent  of  recombination  throughout  the 
expansion  process.  A  brief  description  of  the  Bray  method  for 
determining  the  location  of  departure  from  equilibrium  flow  in  a 
rocket  nozzle  is  given.  The  analysis  is  applied  to  aerozine, 

H2-Oz  and  H2~F2  propellants  and  a  hypothetical  H2~C12  system 
to  obtain  the  point  of  effective  freezing  as  a  function  of  area 
ratio.  The  kinetic  rates  for  the  reactions  occurring  in  the  com¬ 
bustion  process  are  discussed,  and  the  values  for  the  rates  are 
included. 


-  v- 


CONTENTS 


NOMENCLATURE .  ix 

I.  INTRODUCTION .  1 

II.  NOZZLE  PARAMETERS  FOR  BRAY  ANALYSIS .  7 

III.  HYDROGEN-OXYGEN  SYSTEMS  .  9 

IV.  N204-AER0Z INE  SYSTEMS .  15 

V.  HYDROGEN- FLUORINE  SYSTEMS .  19 

VI.  HYDROGEN- CHLORINE  SYSTEMS .  29 

REFERENCES .  31 


FIGURES 


1.  N2C>4-  Aerozine,  Kinetic  and  Equilibrium  Specific 

Impulse . 

2.  Hydrogen- Fluorine  Freeze  Points  . 

3.  Hydrogen- Fluorine,  Kinetic  and  Equilibrium  Specific 

Impulse . 

4.  Temperature  at  Freeze  Point  for  Hydrogen-Fluorine  . 

5.  Recombination  at  Freeze  Point  for  Hydrogen- Fluorine, 

Per  Cent . 


18 

22 

25 

26 

27 


TABLES 


1.  Comparison  of  Frozen  and  Equilibrium  Specific 

Impulses  . 

2.  Gas  Compositions,  Pc  =  1000  lb  in-^  . 

3.  Gas  Composition  for  Hydrogen- Oxygen  Systems, 

P  =  100  lb  in"2  . 

c  .  -2 

4.  N2O4- Aerozine  Systems,  Pc  -  100  lb  in  . 

5.  N2O4- Aerozine  Freeze  Points  and  Performance  Losses 

6.  Hydrogen- Fluorine  Freeze  Points  and  Performance 

Losses  . 


2 

3 

10 

15 

17 

24 


-  Vlll- 


NOMENCLATURE 


A,  B 

,  - 1 

concentrations  of  reactants,  moles  gm 

c* 

-  1 

characteristic  velocity,  ft  sec 

g 

gravity,  32.  17  ft  sec 

/cc  mole  1  sec  for  bimolecular 

reactions 

kb 

backward  rate  constant,  {^2  mol0-2  Se<;-1  for  three-body 

reactions 

kf 

forward  rate  constant 

I 

s 

specific  impulse 

M 

average  gas  molecular  weight,  g  mole 

MR 

mixture  ratio,  weight  of  oxidizer /weight  of  fuel 

N 

number  of  moles  per  reaction,  moles  g 

,,  .  -2 

P 

c 

chamber  pressure,  lb  in 

rk 

kinetic  rate  of  recombination 

r 

s 

recombination  rate  required  for  equilibrium 

rt 

throat  radius,  in. 

R 

gas  constant,  1. 987  cal 

RC 

rate  controlling 

T 

temperature,  °K 

Tf 

temperature  at  freeze  point,  °K 

X 

distance,  ft 

£ 
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NOMENCLATURE  (Continued) 


P  nozzle  half  angle 

«  area  ratio 

«£  area  ratio  of  freezing 

p.  velocity,  ft  sec"' 

p  density,  gm  cc  ' 
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I.  INTRODUCTION 


At  the  high  temperatures  usually  encountered  in  the  combustion  of 
rocket  propellants,  the  product  gases  are  composed  of  stable  molecules, 
partially  dissociated  molecules,  and  free  radicals.  The  calculation  of 
theoretical  propellant  performance  and  gas  composition  by  means  of  shifting 
equilibrium  methods  is  based  on  the  assumption  of  infinitely  rapid  rates  of 
reaction  in  the  combustion  chamber  and  during  the  nozzle  expansion  process. 
Maximum  propellant  performance  is  predicted,  since  at  any  pressure  and 
temperature  the  enthalpy  release  corresponding  to  maximum  recombination 

is  achieved. 

Frozen  propellant  performance  is  based  on  the  assumption  of  zero 
reaction  rates  starting  at  some  point  in  the  rocket  motor.  The  classic 
freeze  point  in  the  older  literature  was  the  chamber  gas  composition.  From 
the  chamber  throughout  the  nozzle  expansion  the  gas  composition  remained 
the  same.  More  recent  approximate  methods  of  calculation  show  that  in 
actual  engines  the  freeze  point  should  be  in  the  vicinity  of  the  throat.  The 
exact  area  ratio  where  freezing  takes  place  varies  with  the  motor  and 
propellant  so  each  case  must  be  interpreted  individually  to  determine  the 
performance  for  the  system. 

A  comparison  of  performance  achieved  from  equilibrium  flow  through¬ 
out  the  nozzle  and  from  freezing  at  different  points  is  shown  in  Table  1  for 
several  typical  propellants. 

If  the  composition  of  the  combustion  products  of  various  propellants 
is  frozen  at  the  throat  or  not  permitted  to  change  as  a  function  of  the  ex¬ 
pansion  process,  the  effect  on  performance  varies  with  the  chemical 
composition  and  the  mixture  ratio  of  the  propellant  (Table  1).  For  example, 
the  Al-H^-O^  propellant  suffers  almost  no  loss  in  performance  when  it  is 
frozen  at  the  throat,  while  OF2-B2H&  shows  a  considerable  loss.  This 
difference  is  obvious  when  the  composition  (Table  2)  at  the  throat  is  examined. 
The  number  of  moles  of  atoms  and  radicals  available  for  recombination 


Table  2.  Gas  Compositions, 


P  = 
c 


1000  lb  in 


-2 


Combustion 

Products 

of2-b2h6 

MR  =  4.  0 

Combustion 

F2-N2H4 

MR  =2.3 

Moles/100  g 

Chamber 

Throat 

Moles/ 1 00  g 

Chamber 

Throat 

H 

1. 03752 

0.  97877 

H 

0.  41974 

0.  34679 

B 

0.  00055 

0.  00031 

N 

0.  00403 

0.  00237 

O 

0. 12013 

0. 09378 

F 

0. 54595 

0.  43794 

F 

0. 12069 

0. 09081 

NH 

0.  00060 

0.  00032 

BH 

0.  00002 

0. 00001 

HF 

3.  2237 

3.  23044 

hbo2 

0.  03885 

0. 04168 

h2 

0.  1  1953 

0.  1021  3 

OH 

0.  13316 

0.  11301 

nh2 

0.  00003 

0.  00001 

HF 

1. 54425 

1. 56241 

n2 

0.  94313 

0. 9441 1 

h2 

0.  67466 

0. 70094 

f2 

0.  00005 

0. 00002 

h2o 

BO 

0.  11543 
0.  11552 

0.  11813 

n  fiQQQi 

H,  A,  F,  NH 

0. 9703 

0.  7874 

BOF 

0.  90003 

0.  95770 

H?/0, 

BF 

bf2 

0. 38122 

0. 33753 

MR  = 

3.  1 

bf3 

0. 00490 

0.  00432 

H 

0. 03350 

0. 01261 

b2o2 

0. 00083 

0. 00073 

O 

0.  00001 

0.  00000 

b2o3 

0. 00044 

0. 00048 

OH 

0.  00414 

0.  00109 

02 

0. 00826 

0. 00635 

h2 

7.  35794 

7.  36683 

H,  B,  O,  F,  OH 

1. 412 

1.  277 

h2o 

o2 

4.  72146 

0.  00000 

4.  72453 

0.  00000 

h2- 

F2 

H,  O,  OH 

0.  03765 

0. 01370 

MR  = 

4.  0 

H 

h2-< 

d2 

0.  11172 

0. 00034 

0. 04129 

0.  00006 

MR  = 

4.  5 

F 

HF 

4. 21019 

4.  21047 

H 

0.  18527 

0. 11961 

h2 

7. 75972 

7.  79476 

O 

0. 00270 

0. 00096 

H,  F 

0.  11206 

0.  04135 

OH 

h2 

0.  08432 

3.  85918 

0. 04506 

3.  86883 

h2-i 

^2 

h2o 

5.  02405 

5.  06669 

MR  = 

7.  0 

o2 

0. 00143 

0. 00052 

H,  O,  OH 

0. 27229 

0. 1656 

H 

0. 68524 

0.  47401 

F 

0. 01881 

0.  00868 

HF 

4. 58645 

4.  59658 

h2 

3. 56458 

3.  66586 

H,  F 

0. 70405 

0.  48269 

Table  2. 

Gas  Compositions,  P 

=  1000  lb  in"2 
c 

(Continued] 

Combustion 
Products 
Moles/ 100  g 

H  0-1 
MR  = 

r2 

11.0 

Combustion 

Products 

H2-02 

MR  =  6.  33 

Chamber 

Throat 

Moles/  1  00  g 

Chamber 

Throat 

H 

F 

HF 

H2  ^ 

H,  F 

1. 08504 

0.  15996 

4. 66476 

1. 25705 

1. 24500 

0. 90394 

0.  09785 

4.  72689 

1. 31656 

1. 00179 

H 

O 

OH 

h2 

h2o 

°2 

H,  O,  OH 

0.  27057 

0.  04155 

0.  37867 

1. 57128 

4.  87270 

0.  05212 

0.  69079 

0. 22593 

0. 02903 

0. 29929 

1. 51476 

4. 99120 

0. 03882 

0.  5542 

AI-H2-O2 

Al-20%  O2/H2  =  2.0 

H 

h2o 

h2 

A1203(c) 

OH 

H,  OH 

0. 01735 

2.  22129 
10. 99398 

0. 37064 

0. 00052 

0.  01787 

0.  00654 

2.  22168 
10. 99919 

0.  37064 

0. 00014 

0.  00668 

Be-H2-Oz 

Be-  2  5%  02/H2  =1.91 

H 

Be 

O 

BeH 

BeOH 

OH 

h2 

Be(OH)2 

h2o 

BeO 

Be202 

Be303 

Be404 

Be£,06 

BeO(c) 

H,  Be,  O,  OH 

0.  16324 

0.  00112 

0.  00001 

0.  00025 

0. 03730 

0.  00076 
12. 38739 

0.  0001  3 

0.  21012 

0.  00001 

0.  00004 

0.  00060 

0.  00000 

0. 00000 
2. 79961 

0.  16513 

0.  16715 

0.  00110 

0. 00001 

0. 00018 

0.  03720 
0. 00072 
12. 38558 

0. 00012 

0. 21008 

0. 00001 

0.  00004 
0. 00082 

0. 00001 

0. 00001 

2. 79909 
0.  16898 

varies  for  the  propellant  and  mixture  ratio.  In  the  Al-H^-C^  propellant 
there  are  only  0.  007  moles  of  H  and  OH  which  would  combine  to  form  H^O 
and  therefore  release  additional  heat  if  the  system  were  permitted  to  proceed 
to  complete  the  recombination.  The  loss  in  energy  by  restricting  the  re¬ 
combination  is  small  and  therefore  causes  only  small  losses  in  performance. 
In  the  OF^-B^H^  propellant  there  are  1.  28  moles  of  H,  B,  O,  F,  and  OH 
which  are  not  permitted  to  recombine,  causing  larger  energy  and  perform¬ 
ance  losses.  Performance  losses  are  proportional  to  the  energy  losses 
caused  by  the  freezing  process. 

A  comparison  of  performance  degradation  among  the  various  mixture 
ratios  of  H^-O^  reveals  the  same  relationship.  At  a  mixture  ratio  of  3.  1, 
only  0.  014  atoms  are  restricted  from  recombining,  while  at  6.  3,  0.  554 
moles  remain  as  free  atoms  or  radicals  after  the  freezing  process. 

Available  specific  impulse  is  determined  by 

iJM,„.h,]1/2  (i) 

s  [ g  c  ej 

where  Ig  is  specific  impulse,  g  is  the  gravitational  constant,  J  is  the 

mechanical  equivalert  of  heat,  Hc  is  the  chamber  enthalpy,  and  Hg  is  the 

exhaust  enthalpy.  In  order  to  obtain  maximum  specific  impulse  it  is  thus 

necessary  to  maximize  H  and  minimize  H  .  The  value  of  H  is  fixed  by 
*  c  e  c 

the  propellant  combination,  but  Hg  is  a  function  of  the  nozzle  expansion  ratio, 
the  kinetics,  and  the  enthalpies  of  the  chemical  reactions  which  occur  during 
expansion.  The  combination  of  Eq.  (1)  with  its  derivative  (Ref.  1)  with 
respect  to  H^  leads  to  Eq.  (2): 


dl  dH 

T~  =  2(H  -  H  T 

s  e  c 
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which  shows  that  the  fractional  loss  in  specific  impulse  is  proportional  to  the 
He  loss  due  to  unrecombined  free  radicals.  In  many  chemical  rockets  this 
specific  impulse  loss  can  be  as  high  as  10  per  cent. 

In  actual  nozzles,  the  residence  time  of  the  combustion  gases  is  of  the 
order  of  a  few  milliseconds,  and  the  idealized  conditions  of  shifting  thermo¬ 
dynamic  equilibrium  may  not  be  achieved.  For  actual  nozzles,  some  re¬ 
combination  reactions  usually  occur,  and  the  achievable  performance  will 
lie  between  the  shifting  and  frozen  equilibrium  values.  The  magnitude  of 
the  nozzle  recombination  performance  loss  will  depend  on  factors  such  as 

Pc*  Po'  €l  and  nozzle  geometry  which  determine  velocity  and  residence 
time  of  the  combustion  gases. 

Because  of  the  complexity  of  a  complete  kinetic  analysis  of  nonequilib¬ 
rium  flow  in  nozzles  and  lack  of  accurate  kinetic  rate  constants,  considerable 
interest  has  been  directed  toward  approximate  methods  of  calculation. 

Bray  (Ref.  2)  in  1959  proposed  that  nozzle  recombination  losses  could  be 
approximated  by  assuming  equilibrium  flow  to  a  selected  point  in  the  nozzle 
followed  by  frozen  flow  for  the  remainder  of  the  nozzle  expansion  process. 

In  a  recombination  reaction  the  principal  mechanism  governing  recombina¬ 
tion  is  a  three-body  reaction,  A+B+M  -  AB  +  M,  in  which  A  and  B  are 
identical  or  dissimilar  free  radicals,  and  M  is  any  third  body.  The  energy 
liberated  in  three-body  recombination  reactions  is  usually  at  least  an  order 
of  magnitude  greater  than  binary  exchange  reactions  of  the  type  A  +  B  —  C  +  D. 
Bray  showed  that  although  there  are  many  possible  three-body  recombination 
reactions  the  number  of  freeze  points1  is  equal  to  the  number  of  free  radicals 
less  the  number  of  independent  binary  exchange  reactions.  Reduction  of  the 
number  of  freeze  points  occurs  because  of  the  rapid  binary  exchange  re¬ 
actions  which  couple  the  free  radical  and  recombination  reactions.  Since  the 
Bray  approximation  method  is  not  a  routine  calculation  but  varies  with  gas 
composition,  the  methodology  is  best  illustrated  by  the  examples  which  follow. 


These  are  the  points  in  the  rocket  nozzle  at  which  equilibrium  gas  composi- 
i°n  cannot:  be  maintained  and  beyond  which  a  frozen  gas  composition  is 
expanded  lsentropically.  * 
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II.  NOZZLE  PARAMETERS  FOR  BRAY  ANALYSIS 


At  the  freeze  point  the  kinetic  rate  and  the  recombination  rate  required 

by  the  nozzle  must  be  of  the  same  order  of  magnitude  to  be  consistent  with 

the  concept  of  the  rapid  transition  from  equilibrium  to  frozen  composition. 

Bray's  detailed  comparison  showed  that  the  ratio  of  the  above  two  rates  at 

the  freeze  point  was  close  to  unity  (Ref  2).  Thus,  he  assumed  that  the  freeze 

point  could  be  approximate d  by  equating  the  kinetic  rate,  r  ,  to  that  required 

bv  the  nozzle,  r  . 

s 

The  reaction  rate  required  to  maintain  shifting  equilibrium  in  the 
nozzle  is  given  by  Eq.  (3): 


2gP  tan  (3 

i-  -  c  1 

fdNi 

1780.  7P  /.  _x 

I  -  ,  cftan  P] 

/dN\ 

rs  r  ' 
s  C*r  .\T« 

Ide  J 

£*rt  \  n r e  / 

\d£  ) 

(3) 


For  the  derivation  of  Eq.  (3)  it  is  assumed  that  the  reaction  rate 
required  to  maintain  shifting  equilibrium  in  the  nozzle  may  be  defined  by 
the  following  equation 


where  p  is  density,  p  is  gas  velocity,  and  x  is  linear  nozzle  distance.  From 
geometrical  considerations  it  can  be  shown  that 

de  _  Z\Ti  tan  (3 
d^  -  rt 


and  from  the  theory  of  rocket  propulsion  that 

gP 


=7c? 
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In  the  preceding  equations  (3  is  the  nozzle  half  angle,  is  the  nozzle  throat 
radius,  g  is  the  gravitational  constant,  Pc  is  chamber  pressure,  and  C*  is 
the  characteristic  velocity.  Substitution  of  the  relationships  for  dc/dx  and 
pp  into  the  equation  for  r^  leads  to  Eq.  (3). 

Solution  of  Eq.  (3)  requires  evaluation  of  dN/de  from  the  shifting 
equilibrium  nozzle  compositions  and  tan  p/\/e  by  measuring  the  slope  at 
various  e  stations  in  the  nozzle.  For  a  conical  nozzle  p  is  constant,  except 
for  the  curved  section  near  the  throat.  Substitution  of  concentrations  of  the 
reactants  at  various  values  of  t  permits  the  calculation  of  r^,  the  kinetic 
rate  of  combination.  For  a  three-body  reaction  r^  =  k^ABM  when  con¬ 
centrations  are  expressed  in  moles  per  cc  or  k^ABp^/M  when  concentra¬ 
tions  are  expressed  in  moles  per  gram,  and  M  is  the  average  molecular 
weight  of  the  exhaust  gases.  The  third-body  concentration  is  taken  to  be  the 
total  gas  composition.  The  value  of  e  at  which  r  =  r,  is  by  definition  the 
freeze  point.  The  value  of  the  kinetic  specific  impulse  is  determined  by 
employing  shifting  equilibrium  composition  to  the  freeze  point  and  freezing 
the  gas  composition  for  the  remainder  of  the  nozzle  expansion  process 
(Ref.  3). 
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III.  HYDROGEN- OXYGEN  SYSTEMS 


The  first  step  for  any  Bray  analysis  is  the  determination  of  the 
important  chemical  species  formed  from  the  combustion  process  and  the 
pertinent  chemical  reactions  occurring  between  the  species.  When  possible, 
the  equations  are  written  in  an  exothermic  direction,  and  the  products  of 
the  reactions  are  those  whose  concentrations  increase  during  the  expansion 
process. 

The  changes  in  equilibrium  gas  composition  which  occur  as  a  function 
of  mixture  ratio  are  illustrated  in  Table  3.  For  hydrogen- oxygen  systems 
with  mixture  ratios  as  low  as  1.  6  and  2.  4,  such  small  amounts  of  free 
radicals  are  present  that  recombination  reactions  are  not  significant.  The 
composition  is  almost  frozen  in  the  chamber. 

At  the  higher  mixture  ratios,  there  are  sufficient  free  radicals  present 
to  influence  the  energy  released  from  possible  recombination  reactions.  It 
now  becomes  practical  to  estimate  a  Bray  freeze  point  by  choosing  the 
kinetic  steps  which  apply.  Since  OH,  H  and  O  are  free  to  form  H^,  H^O  and 
02,  reactions  involving  these  radicals  are  considered.  Possible  steps 
might  consist  of  the  following  (Ref.  4). 
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Since  there  are  six  species  containing  two  elements,  four  reactions  con¬ 
taining  all  six  species  are  needed  to  describe  or  couple  the  changes  in 
concentration.  Both  H2  and  HzO  increase  throughout  the  expansion  process, 
while  H,  OH,  O,  and  02  decrease.  Wherever  possible  it  is  desirable  to 
write  the  equation  in  the  direction  of  increasing  H2  and  H20  to  describe  the 
chemical  process  actually  occurring  and  also  to  maintain  a  positive  concen¬ 
tration  change  in  the  final  algebraic  expression. 

Assuming  unit  concentrations,  the  possible  reaction  can  be  inspected 
for  their  relative  rates  and  for  their  energy  changes.  The  bimolecular 
exchange  reactions,  ©  through  ©  ,  explain  the  interdependence  of  the 

chemical  species.  They  are  faster  than  the  three  body  reactions,  © 

through  ©,  and  therefore  are  assumed  to  be  more  nearly  m  equilibrium 
and  are  not  rate- controlling  phenomena.  The  three-body  reactions  are  very 
exothermic  and  represent  the  recombination  occurring  in  the  nozzle  expan¬ 
sion.  Since  the  three-body  reactions  release  the  most  energy  and  are  rate 
controlling,  the  pertinent  ones  should  be  selected.  For  the  hydrogen- oxygen 
case  the  water  and  hydrogen  recombination  steps  should  be  considered. 

Both  O  and  OH  decrease  in  concentration  with  time,  so  their  terminal 
three-body  reactions  are  discounted.  Two  other  equations  from  ©  through 
©  should  be  selected  to  adequately  account  for  all  the  species.  Reaction 
©is  not  kinetic  ally  feasible  process  because  of  the  high  activation  energy 
needed  to  break  the  H.,  and  O.,  bonds.  The  remaining  reactions  have  similar 
rates,  assuming  unit  concentrations,  and  a  choice  among  them  becomes 
mostly  an  algebraic  one.  The  four  reactions  chosen  were 

©  H+H+M  —  H2+M 
©  H  +  OH  +  M  -  HzO  +  M 
©  C>2  +  H  -*■  OH  +  O 
©  HzO  +  O  -  2  OH 


-11- 


In  describing  each  reaction  in  terms  of  the  formation  or  disappearance 
of  one  specie  the  following  expressions  are  obtained: 


dH2° 

de 


+  r; 


where 


r 


1 


dH 

de 


2r2- 


r6 


and  where 


dH. 


d°2 

dc 


Substitution  into  the  above  expressions  yields 


dH  2dH2  +  dG2 


de  de 


de 


d(-H20-  H-  2H2  +  02)  _dN 
de  de 


The  choice  between  the  two  three-body  reactions,  (T)  and  (2)  ,  to 
determine  the  rate- controlling  step  to  be  used  in  the  Bray  analysis  is  some¬ 
what  more  difficult.  One  approach  is  the  substitution  of  actual  concentrations 

2 

of  the  participating  species  in  the  equation  r^  =  k^ABp  for  bimolecular  reactions 
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or  r,  =  kf 
k  f 


ABp3/M. 


The  density  is  determined  from  the  gas  law, 


PM 
p  =  RT 


where 

p  =  density,  g  cc  ^ 

P  =  pressure,  atm 

M  =  average  gas  molecular  weight,  g  mole 

T  =  temperature,  °K 

R  =  gas  constant,  82.  057  atm  cc  deg  1  mole  1 

The  evaluation  of  r,  for  the  throat  conditions  yielded  a  value  for  reaction 
?  k  _  3 

(2)  of  1 . 4  X  10"  ,  while  r^  for  reaction  QT)  was  2.8X10  .  Reaction 

^1)  appears  to  be  rate  controlling  by  this  criterion,  but  since  the  concen¬ 
tration  influences  the  final  result,  a  concentration  of  little  importance  to  the 
total  chemical  process  could  be  construed  to  control  the  entire  process. 

For  this  reason  the  changes  in  H2  and  HzO  concentration  in  a  equilibrium 
process  might  be  considered.  At  a  mixture  ratio  of  4.  5  the  AHjO  =  0.  069 
and  the  AH^  =  0.013  moles/100  g  between  the  chamber  and  the  throat.  This 
relative  difference  continues  beyond  the  throat;  therefore,  a  change  in  the 
formation  would  influence  the  performance  more  than  a  change  in  H2 
formation.  The  enthalpy  changes  in  the  two  reactions  are  almost  equal.  As 
a  first  approximation  H  +  OH  —  H20  looks  like  a  rate- controlling  step. 

In  fact,  H+H  +  M— H2  +  M  might  be  the  rate- controlling  step.  Since 
the  two  reactions  have  similar  rates,  two  freeze  points  could  be  determined, 
using  each  of  the  rate- controlling  steps  separately. 

Evaluation  of  the  freeze  point  for  the  same  nozzle  for  a  hydrogen- 
oxygen  propellant  at  a  mixture  ratio  of  4.  5  and  chamber  pressure,  Pc»  of 
100  gave  only  slightly  different  freeze  points  by  varying  the  rate- controlling 
step.  For  the  H+H  +  M  —  H2  +  M  step  the  freeze  point  occurred  at  an  area 
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ratio  of  1. 26,  representing  a  1.  3  per  cent  loss  in  vacuum  specific  impulse 
at  an  area  ratio,  e ,  of  50.  When  H  +  OH  +  M  -  HzO  +  M  was  considered  as 
the  rate- controlling  step,  the  freeze  point  was  at  an  area  ratio  of  1. 58,  which 
gave  an  impulse  2.  5  sec  higher  than  the  area  ratio  of  freezing,  e^,  of  1. 26. 
The  temperatures  at  the  freeze  points  were  2568  and  2424°K,  respectively. 

The  kinetics  developed  in  the  following  section  for  N204~  aerozine 
were  also  applied  to  the  same  H2-C>2  propellant  and  yielded  a  freeze  point 
of  1.  39,  an  intermediate  value  between  the  two  points  found  using  the  kinetics 
developed  in  this  section.  For  H2~02  the  variations  in  the  selected  kinetics 
produced  only  small  changes  in  the  freeze  point.  Since  the  smallest  value, 

1.  26,  represents  the  largest  impulse  loss,  it  is  the  most  conservative 
estimate  of  the  freeze  point. 


IV.  N_.0  .-AEROZINE  SYSTEMS 
2  4 


The  composition  of  an  N^-  aerozine  propellant  at  a  mixture  ratio  of 
2  is  shown  in  Table  4.  The  analysis  which  follows  is  valid  for  mixture 
ratios  between  1.  6  and  2.  3. 


Table  4.  N2C>4- Aerozine  Systems,  Pc  -  100  lb  in 


MR  = 

l.  0 

Composition,  Moles/100  g 

Chamber 

Throat 

H 

0. 0935 

0.  07376 

o 

0. 03898 

0.  02866 

OH 

0. 20515 

0.  16779 

H,0 

1. 70138 

1.  7615 

2 

H._ 

0. 29036 

0.  25881 

2 

CO 

0. 29099 

0.  26668 

CO, 

0. 23034 

0.  25466 

NO 

0. 0504 

0. 03929 

N_ 

1.49253 

1. 49810 

2 

°2 

0. 09179 

0.  07897 

The  kinetic  steps  given  for  the  H.,-0-,  system  apply,  and  in  addition 
the  reaction  CO  +  OH  -*•  CO.,  +  H  should  be  considered.  The  k  for  this  re¬ 
action  is  expt-104  cal/ RT)  cc  mole'1  sec  and  the  AH^g  -  -24.88  kcal. 
The  reaction ‘is  fast  compared  to  the  three-body  reactions  and  couples  with 
3  +  H  -  H  +  OH  reaction,  giving  the  water  gas  reaction  with  a  net 
of  -9.  8  kcal.  The  total  energy  gain  is  small.  Nitrogen  interacts 


the  H  O  +  H  —  H  +  OH  reaction,  giving  the  water  gas  reaction  with  a  net 

2  ^  -j.  _ 

AH 

with’S.e"  oxygen  and  hydrogen  specie,  to  such  a  small  extent  that  it  is  con¬ 
sidered  inert.  The  only  remaining  specie,  to  include  are  the  hydrogen- 
oxygen  combinations  which  are  the  same  a,  the  product,  in  the  hydrogen- 

oxygen  case. 
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V 


The  greatest  difference  between  the  combustion  products  of  the  two 
propellants  is  the  change  in  H2  concentration.  For  H2-02>  H2  increases 
down  the  nozzle,  while  it  decreases  from  the  chamber  to  an  «  of  8  for 
N204-aerozine.  Since  this  is  the  area  where  freezing  is  expected  to  occur, 
the  H  +  H  +  M-*H2+M  reaction  no  longer  applies.  This  leaves 
H  +  OH  +  M  —  H20  +  M  as  the  rate- controlling  step;  therefore,  only  one 
freeze  point  needs  to  be  considered.  The  final  set  of  equations  chosen  for 
the  analysis  were 


@  H+OH+M-  H20  +  M 
©  02  +  H  +  OH  +  O 
©  O  +  H20  +  2  OH 
©  H2  +  OH  +  H  +  HzO 

The  expression  of  the  rates  as  concentration  changes  results  in  the 
following  relationships: 


rl  = 


dH2° 


r  3 


d  e 

dO 


+  r. 


de  +  r2 


r4  =  - 


dH2 

de 

dO_ 


r2 
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Substitution  into  the  above  expressions  yields 


r3  d«  d« 

d(HzO  -  O  -  Oz  +  Hz)  _  dN 
rl  “  de  dt 

A  Bray  analysis  was  calculated  with  these  kinetics  to  determine  the 
freeze  points  for  four  mixture  ratios.  Table  5  shows  the  freeze  points, 
kinetic  impulses,  and  performance  losses  (per  cent)  when  compared  to  equi¬ 
librium  impulses.  The  same  results  are  shown  graphically  in  Fig.  1. 


Table  5.  H2O4-  Aerozine  Freeze  Points 

and  Performance  Losses 


RC  Step 

P 

C 

MR 

c £»  Trans 
Nozzle 

* 

H 

Is(vac) 

€  =  10 

Loss 

from  equil,  % 

Is(vac) 

«  =  50 

Loss 

from  equil,  % 

H  +  OH  +  M  - 

100 

1.6 

1.  28 

2627 

305.  2 

1.2 

329-  5  • 

1.7 

HzO  +  M 

1 .  8 

1.  43 

2674 

305.  1 

1 . 6 

330.  4 

2. 4 

2.  0 

1.  48 

2701 

302.  3 

2.  2 

327.  8 

4. 0 

2.  3 

1 . 40 

2711 

295.  1 

2.  5 

319.  9 

4.  8 

- - - 
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zine,  Kinetic  and  Equilibrium  Specific  Impulse 


V.  HYDRO  GEN- FLUORINE  SYSTEMS 


The  products  formed  from  a  hydrogen- fluorine  propellant  are  H^,  HF, 
H,  and  F.  The  product  F^  is  present  only  at  extremely  high  mixture  ratios 
and  then  in  small  concentrations  which  are  unimportant  for  Bray  analysis. 
The  reactions,  their  rate  constants,  and  AH's  selected  for  consideration  are 
as  follows: 


The  reaction  —  2  HF  is  not  considered  a  likely  reaction 

kinetically.  The  product  HF  is  formed  by  one  of  the  above  mechanisms. 

The  three-body  reactions,  (T)  and  (2),  represent  the  possible  rate¬ 
controlling  steps,  since  (5)  is  a  faster  bimolecular  reaction  releasing  less 
energy  and  is  more  nearly  in  equilibrium  than  (7)  or  (Z)  . 

The  product  HF  was  chosen  as  the  third  body  since  it  is  present  in 
larger  quantities  than  any  other  gas  and  is  the  largest  contributor  to  the 
average  molecular  weight,  which  was  used  in  the  actual  calculations  as  the 
third-body  concentration. 


A  point  to  consider  in  choosing  the  rate- controlling  step  is  whether 
or  not  it  represents  the  chemical  phenomena  occurring  in  the  nozzle.  For  a 
hydrogen-fluorine  propellant  at  a  chamber  pressure  of  150  psi  and  a  mixture 
ratio  of  12,  there  are  1.  27  moles  H,  0.  17  moles  F,  4.  68  moles  HF,  and 
0.84  moles  in  100  grams  for  the  throat  conditions  of  a  shifting  equilibrium 
calculation.  Since  chamber  pressure  changes  composition  very  little,  the 
analysis  which  follows  is  applicable  to  the  other  pressures  considered. 

At  a  mixture  ratio  of  12  or  18  the  H  +  F  +  M  -  HF  +  Jvf.  is  the  pre¬ 
dominant  reaction.  Upon  substitution  of  the  concentrations  of  the  combustion 
products  in  an  equilibrium  calculation,  the  HF  reaction  is  slower  than  the 
^2  r ecombination  and  therefore  was  chosen  as  the  rate- controlling  step. 

Choosing  two  independent  rate  equations  and  expressing  them  in  terms 
of  the  species  and  HF ,  leads  us  to  the  following  equations: 


_  d(HF) 
2  dt 


d(H  ) 

and  r  3  =  '  -37“ 


Upon  substitution 


_  d(HF)  v  d(H2)  _  d<HF  +  H2> 

2  d  €  d  €  d  c 

This  final  expression  of  r^  is  used  for  the  nozzle  recombination 
dN/df . 

This  particular  set  of  equations  is  valid  for  mixture  ratios  which 
yield  a  large  amount  of  F  and  H  which  may  recombine  into  HF  and  How¬ 

ever  for  lower  mixture  ratios,  4  to  8,  which  are  deficient  in  free  fluorine, 
it  is  a  dubious  treatment.  The  chemical  process  occurring  is  primarily 
the  recombination  of  hydrogen  which  is  expressed  by  and  should  be  used 
as  the  rate- controlling  step. 
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For  comparative  purposes  the  freeze  point  was  determined  for  mixture 
ratios  5.  2  and  8.  0,  using  k1  and  k2  separately  as  the  rate- controlling  steps. 
The  freeze  points  and  corresponding  specific  impulses  are  shown  in  Table  6. 

At  a  mixture  ratio  of  8,  for  example,  the  use  of  rates  k1  and  k2  yields 
freeze  points  of  1.  50  and  3.  85  for  the  transtage  nozzle,  respectively,  cor¬ 
responding  to  I  losses  of  2.  0  and  0.  5  per  cent  at  an  area  ratio  of  50.  The 
absolute  difference  in  specific  impulse  is  8  sec.  The  rate- controlling  step 
selected  can  influence  the  freeze  point  significantly  and  therefore  should  be 
selected  carefully.  In  general,  the  larger  the  freezing  area  ratio,  the 
smaller  the  performance  loss. 

Since  the  rate  of  the  H  +  F  +  M  -  HF  +  M  was  estimated,  the  effect  of 
a  different  value  was  investigated.  The  constant  was  decreased  one  decade 
to  1  X  10^T-^'  At  a  mixture  rate  of  12  the  freeze  point  dropped  from 
2.  99  to  1.  79,  using  the  new  rate.  The  performance  varied  by  7.  6  sec  at  an 
area  ratio  of  10.  Similar  losses  occurred  for  the  other  mixture  ratios  and 
are  noted  in  the  tabulation.  Until  the  reaction  rate  constants  are  determined 
experimentally  it  is  difficult  to  find  an  exact  freeze  point  from  the  Bray 
method.  Various  alternatives  must  be  considered  to  obtain  minimum  and 
maximum  values. 

A  comparison  of  performance  losses  due  to  kinetic  processes  for 
different  chamber  pressures  shows  that  the  losses  are  greater  at  the  lower 
pressures  since  freezing  occurs  at  smaller  area  ratios.  For  a  H2~F2 
propellant  at  a  mixture  ratio  of  12  and  e  of  10  there  is  only  0.  95  per  cent 
performance  loss  for  a  transtage  nozzle  at  a  pressure  of  300  while  the  loss 
is  2.  7  per  cent  at  100  psi. 

The  kinetic  freeze  point  occurs  at  higher  area  ratios  for  higher 
mixture  ratios  for  the  H2-F2  series  (Fig.  2).  This  seems  logical  since  the 
total  number  of  uncombined  H  and  F  atoms  in  the  chamber  increases  with 
mixture  ratio.  At  a  mixture  ratio  of  5.  2  there  are  0.  48  moles  of  atoms  and 
at  18  there  are  1  .  78  moles.  Essentially,  the  HF  and  H2  recombination 

processes  are  already  more  complete  in  the  chamber  at  low  mixture  ratios; 
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therefore,  the  freeze  point  occurs  nearer  the  chamber.  This  does  not 
prove  that  performance  losses  are  the  least  at  low  mixture  ratios,  since 
the  freeze  point  moves  downstream  and  allows  the  energy  producing  re¬ 
actions  to  occur  at  the  higher  mixture  ratios.  Specific  impulses  and  per¬ 
formance  losses  are  shown  in  Table  6  and  Fig.  3. 

Since  a  complete  Bray  analysis  for  a  series  of  propellants  is  a  tedious 
calculation,  we  attempted  to  find  a  simpler,  empirical  correlation.  Assuming 
the  same  rate- controlling  step  for  each  mixture  ratio,  we  found  that  the 
temperature  at  each  freeze  point  was  similar.  Although  the  chamber  tem¬ 
perature  for  mixture  ratios  5.  2  through  18.  0  varied  by  1400  deg,  the 
freezing  temperatures  for  a  given  nozzle  varied  only  300  deg.  In  Fig.  4 
the  temperature  versus  area  ratio  for  each  mixture  ratio  is  plotted,  and  the 
freeze  point  is  noted  for  each  of  two  nozzles.  The  freeze  points  for  mixture 
ratios  5.  2  and  8.  0  are  those  which  occur  when  the  hydrogen  recombination 
step  is  considered  the  rate- controlling  step.  The  temperature  at  the  freeze 
points  obtained  from  different  rate- controlling  processes  does  not  fall  within 
the  same  range.  The  various  kinetic  processes  used  in  the  analysis  not  only 
produce  quite  different  freeze  points  but  also  widely  varied  temperatures. 

As  a  first  approximation,  the  freeze  point  for  a  chemical  system 
governed  by  the  same  kinetic  process  as  in  a  previously  analyzed  system, 
should  fall  within  the  same  temperature  range. 

Another  empirical  correlation  was  attempted  by  comparing  the  per¬ 
centage  of  atoms  in  the  chamber  which  were  recombined  at  the  freeze  point, 
determined  from  H+  F  +  M-HF  +  M  reaction,  for  the  various  mixture 
ratios.  Figure  5  shows  the  results:  the  degree  of  recombination  decreases 
for  higher  mixture  ratios  or  decreases  for  the  more  highly  dissociated 
systems . 
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Recombination  at  Freeze  Point  for 
Hydrogen-Fluorine,  Per  Cent 


VI.  HYDROGEN- CHLORINE  SYSTEMS 


In  order  to  compare  the  freeze  point  as  determined  by  the  Bray 
analysis  for  propellants  containing  different  halogens,  an  arbitrary  H^-Cl-, 
propellant  was  selected.  The  equilibrium  performance  was  calculated  for  a 
mixture  ratio,  22,  which  contained  the  same  mole  ratio  of  C1-,  to  H-,  as  an 

H2"F2  ProPellant  at  a  rmxture  ratio  of  12. 

The  chlorine  propellant  contained  only  0.00723  moles  of  atoms  in  the 
chamber  and  0.00148  moles  at  the  throat  which  would  be  available  for  a 
recombination  process.  The  comparable  fluorine  propellant  contained 
1.80  moles  of  atoms  in  the  chamber.  For  all  practical  purposes,  an  inspec¬ 
tion  of  composition  showed  that  the  chlorine  propellant  selected  was  frozen 
in  the  chamber. 

Of  course,  other  chlorine- containing  propellants,  such  as  a  solid 
propellant  with  NH4C104>  could  have  been  selected,  but  the  products  of 
combustion  would  be  too  numerous  for  a  simple  Bray  approach  as  we  were 
applying.  The  following  reaction  rates  were  selected  or  estimated  to  attempt 
a  Bray  analysis. 


Reaction 

Rate  Constants 

AH298^kcal  mole 

© 

H  +  H  +  M-H2  +  M 

k1  =  i  xio19t"0-5 

-104. 2 

© 

H+HC1 -H2 +C1 

k^.-^xlo'V05^0™-1 

-1.21 

© 

H2  +  Cl -HC1 +  H 

k3  =  2.82xlOUT-0-5  -500R°Tcal 

+  1. 21 

© 

H+Cl  +  M-*HC1  +  M 

k4=ix  io19t'°‘  5 

-  103. 0 
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Reactions  (2)  and  (3)  were  obtained  from  Benson  (Ref.  5).  Reactions  CD 
and  were  estimated  from  the  fluorine  data,  assuming  that  HC1  as  the 
third  body  would  affect  the  reaction  in  the  same  way  as  HF,  and  that  the 
H  +  Cl  recombination  is  the  same  as  the  H  +  F  recombination. 

It  was  found,  as  had  been  anticipated,  that  freezing  occurred  before 
the  throat  when  the  Bray  technique  was  applied. 
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